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Summary 




A mathematical model of one experiment ideally 
mounted to the experiment support tube is considered. A 
set of compatibility equations is developed and expanded to 
determine expressions for interaction loads and rotational 
deflections. Two digital computer programs are developed 
to undertake a host of thermal inputs for various experi- 


ment sizes and mounting locations. 
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Symbols 


{S T } temperature induced displacement and rotation matrix for the experiment 

{s T } temperature induced displacement and rotation matrix for the hexagonal support tube 

[f] influence coefficient matrix for loads applied at the neutral axis of the experiment 

[ f] influence coefficient matrix for loads applied at the neutral axis of the hexagonal support 

tube 

{x} matrix of loads acting off the experiment neutral axis and referenced to the neutral axis 
of the experiment 

{x} matrix of loads acting off the hexagonal support tube neutral axis and referenced to the 
neutral axis of the hexagonal support tube 

rotational deflection of point P with respect to point P' , radians (Figure 3) 
rotational deflection of point P" with respect to point P' , radians (Figure 3) 

8 0 p'/p rotational deflection of point P" with respect to point P, radians (Figure 3) 

L distance between forward experiment mount and assumed cantilevered point, in 
V length of hexagonal support tube subjected to free bending, in 

The following symbols, where capped with a bar (-), refer to the hexagonal support tube, 
a coefficient of linear thermal expansion of experiment housing, in/in/°F 
E modulus of elasticity of experiment housing, psi 
I cross sectional moment of inertia of experiment housing, in 4 
T uniform temperature of experiment referenced to a datum temperature, °F 

AT transverse linear thermal gradient applied across experiment, °F 

A cross sectional area of experiment housing, in 2 
h depth of experiment housing, in 

X H horizontal reaction induced at forward experiment mount, lb 
X v vertical reaction induced at forward experiment mount, lb 
X M moment induced at forward experiment mount, in-lb 
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A PRELIMINARY THERMAL DISTORTION ANALYSIS 

OF THE 

PROPOSED AOSO EXPERIMENT SUPPORT STRUCTURE 

by 

Donald E. Gray 


INTRODUCTION 

One of the primary objectives of the Advanced Orbiting Solar Observatory (AOSO) structural 
subsystem is to satisfy and maintain the stringent optical alignment requirements established be- 
tween sun sensors and certain of the contemplated experiments. The spacecraft substructure 
pertinent to this objective is the central experiment support tube. 

The proposed experiment support tube is an aluminum alloy hexagonal cylinder composed of 
six joined honeycomb sandwich panels. The hexagonal tube is joined to the external observatory 
structure through three radial pylon supports which are attached to the tube at its midspan. Pro- 
posed experiment mounting to the support tube is through two yoke end mounts which are joined by 
four diagonal struts. 

The optical misalignment induced between sun sensors and individual experiments is a func- 
tion of both geometric and thermal variables. For an assumed geometric configuration and zero 
"g" environment, the optical misalignments can be attributed primarily to the temperature gradients 
over, and temperature levels of, both the support tube and the experiments. 


DISCUSSION 

It is of first interest to obtain some estimate of the interaction loads induced at the experiment 
mounts by experiment interaction with the support tube. As shown in Figure 1, a single experiment 
is considered mounted to the top of the support tube and left of its midspan. The thermal gradients 
applied across the experiment and tube are assumed linear and constant over the length of the 
structure between mount locations. The Euler- Bernoulli theory is assumed and the experiment 
and tube are considered cantilevered at the aft (right) mount. It is further assumed that no rela- 
tive rotations or displacements occur at the mounts. 

Considering the deflections of points P andP (see Figure 1) and applying the above assumptions, 
a set of compatibility equations can be developed in terms of interaction loads at the forward 
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experiment mount. The sign convention used for deflections, loads, and temperature gradients is 
depicted in Figure 2. Also shown are the loads at points P and p. In matrix form, the com- 
patibility equation is 


i 


{s T } + [f] {x} = {S r } + [f] {x} (1) 

Since X H = -X H , X v = -X v , and X M = -X M , Equation (1) can be expanded to 
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The individual terms of Equation (2) are developed in Appendix A. 

The influence coefficient matrices are asymmetrical due to the non- zero terms f 12 and f 
Both terms must be considered since points P and P are at substantial distances from their re- 
spective neutral axis. 

Adjusting the influence coefficient matrices in Equation (2), the load column matrix can be 
rewritten in terms of the three primary reaction loads X H , x v , and X M . Thus 
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Assuming both the support tube and the experiment housing to be constructed primarily of the 
* same material, E - E and a a a . The load matrix can now be isolated by grouping terms and 
factoring. Thus 



A digital computer program written in Fortran n language for the solution of Equation (4) is 
listed in Appendix B as Program I. Using this program, interaction loads can be calculated for 
any number of cases. Loads for eight example cases were calculated and are listed in section A 
of the results. 

Since optical alignments must be held within required limits, it is also of interest to obtain 
some estimate as to the amount of misalignment induced during flight conditions. The optical 
misalignment can be considered to be the relative angular displacement between two pre- aligned 
points. As shown in Figure 3, a single experiment is considered mounted symmetrically about the 
support tube midspan. The thermal gradients applied across the experiment and tube are assumed 
linear and constant over the length of both structures. Again the Euler-Bernoulli theory is assumed, 
but the experiment and tube are considered cantilevered at the tube midspan. It is also assumed 
that no relative rotations or displacements occur at the mounts. 

Referring again to Figure 3 and selecting P' as a reference, the rotational deflection of Pis 
simply the rotational terms for the hexagonal tube in Equation (2). Thus 


r 


s ev 


-aATL 


[0 -L 2 / 2EI L/Elk 


X„ + X„7 


(5) 


Expanding Equation (5) we get: 
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a p 


-aATL L 2 X v LX m Lh X H 
h + 2EI " fl " 2EI 


or 


& 9 P 


-aATL L /LXv 

h + El \ 2 - X m - 2 ; 


(6) 
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Since there are assumed no relative rotations or displacements at the mounts, Equation (6) also 
represents the rotational deflection of point p with respect to either P' or P\ 

Next consider the rotational deflection of point P" at the left end of the support tube. If point 
P is selected as a reference, the rotational deflection of P"is due only to the temperature gradient 
across the tube. Therefore, if L is replaced by L 7 in the temperature induced rotational displace- 
ment term for the hexagonal tube in Equation (2), we get: 

S 6>p"/p “ -oATL ' / h (7) 

The total rotational deflection of P" with respect to the midspan, P' or P', is equal to the 
algebraic sum of the relative deflections. Thus 


h ep" ~ 8 9p + h ep"/p' 


( 8 ) 


Substituting for S^p and and collecting terms 
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dp" 


-aAT , L / LX v hX H \ 

h + L ^ + El \ 2 Xm 2 / 


(9) 


Since the optical axis of the assumed experiment package has not been defined, the solutions 
of all three Equations (6), (7), and (9) may be of interest. A digital computer program written in 
Fortran II language for the solution of these three equations is listed in Appendix B as Program II. 
The program also yields the interaction loads introduced at points P and P, i.e., the solution of 
Equation (4). 


RESULTS 

A. Interaction Loads Introduced at the Forward Experiment Mount 

Referring to Figure 1, the interaction loads introduced at points P and P were obtained through 
the solution of Equation (4) for the following eight cases: 


Geometrical Properties Utilized 


Hexagonal Tube 

Experiment Housing (Assumed) 

I = 28.03 in 4 * 

I = T/2 and in 4 

A = 1.15 in 2 

A = 0.50 in 2 

h = 14 in 

h = 7 in 

E = 10 7 psi 

E = 10 7 psi 

a = 12.8 X 10“ 6 in/in/°F 

a = 12.8 x 10” 6 in/in/°F 

*See Appendix C for calculations. 
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Temperature Distributions 


(a) 

AT = 0°F 
AT = 1°F 
T = 0°F 

T = 1°F 


(b) 


(c) 


AT = 1°F 
AT = 1°F 
T = 0°F 
T = 0°F 


AT = 0°F 
AT = 3°F 
T = 0°F 
T = 1°F 


Example Cases with Results 


(d) 

AT = 1°F 
AT = 3°F 
T = 0°F 
T = 1°F 


Case 1 - Assume I = l/2, L = 5.75 in, and temperature distribution (a). 
X H = 11.66 lb 
X v = 0. lb 
X„ = -85.45 in-lb 


Case 2 - Assume I = I /2, L = 28.5 in, and temperature distribution (a). 
X H = 11.66 lb 
X v = 0. lb 
X„ = -85.45 in-lb 

Case 3 - Assume I = T/2, L = 28.5 in, and temperature distribution (b). 
X H = -23.311b 
X v = 0. lb 
X* = 85.46 in- lb 


Case 4 - Assume I = °° , L = 5.75 in, and temperature distribution (a). 
X H = 44.61 lb 
X v = 0. lb 
X M = -568.52 in-lb 

Case 5 - Assume I = ®, L = 28.5 in, and temperature distribution (a). 
X H = 44.61 lb 
X v = 0. lb 
X„ = -568.52 in-lb 


Case 6 - Assume I = 00 , L = 28.5 in, and temperature distribution (b). 

X H = -66.91 lb 
\ = 0. lb 
\ = 724.64 in-lb 

Case 7 - Assume I = 1/2, L = 28.5 in, and temperature distribution (c). 
X H = -11.65 lb 
X v = 0. lb 
X„ = -256. 33 in-lb 

Case 8 - Assume I = °°, L = 28.5 in, and temperature distribution (d). 

X H = -22.301b. 

X v = 0. lb 

X„ = -100.15 in- lb 
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It is of interest to note that the vertical force, x v , is equal to zero for all eight cases. For 
the assumed configuration, the vertical force, X v , is always zero. This relation is shown in Ap- 
pendix A. 

B. Induced Rotational Deflections 

Equations (6), (7), and (9) were solved for the configuration depicted in Figure 3. The first 
run utilized the f Blowing input values assuming aluminum alloy for both the experiment housing 
and the hexagonal tube: I = 28.03 in 4 , I = T/8, A = 1.15 in 2 , A = 0.50 in 2 , h = 14 in, h = 7 in, 

E = 10 7 psi, a = 12.8 x 10' 6 in/in/°F, L = 10.51 in, and L' = 16.17 in. A total of 144 different 
cases of temperature-gradient combinations were considered. In many instances, the same de- 
flection resulted for different temperature- gradient combinations. Therefore, cases yielding 
like deflections were grouped together where each group number identifies the first case in that 
particular group. These groups are listed in Table I along with the corresponding temperature- 
gradient values for each case. 

The results are presented as sets of parametric curves for the various groups in Figures 4, 

5, and 6 with deflections plotted as a function of the temperature gradient across the hexagonal 
tube. Figure 4 shows the variation of h 6 p with AT. This deflection is represented in Program II 
by the variable name, DEFHEX (see Appendix B). Since the curve slopes are equal and constant 
for all groups, the effect of changing the temperature- gradient combinations is to shift the curve 
along the ordinate to a new deflection for a given AT reference. This shows that for a given de- 
flection there exists a range of possible AT’s. It is of interest to note that zero deflections are 
theoretically possible for AT’s in the range of ±7 seconds of arc. The apparent subgroups of three 
are for cases containing similar T and T parameters. The variable within these subgroups is the 
less significant at parameter. The single curve in Figure 5 indicates the variation of S e - ypwith 
AT for all groups since, from Equation (7), the free bending of the beam is dependent only on 
the thermal variable, AT, holding a, L', and h constant. This deflection is denoted by the variable 
name DEFREL in Program II. The slope of this function is shown to be -3.051 arc sec/°F. The 
variation of S 0 p„with AT is shown in Figure 6. This deflection is the algebraic sum of DEFHEX 
and DEFREL and has the variable name DEFEND in Program II. As indicated, zero deflections 
are also theoretically possible for AT’s in the range of ±7 seconds of arc. The slope of these 
curves is shown to be -3.845 arc sec/°F. 

Several additional runs were executed varying such parameters as the type of structural ma- 
terial, length of support tube, distance between experiment mounts, and experiment housing cross 
sectional moment of inertia. In addition, runs were executed for zero and exceptionally large tem- 
perature gradients across the experiment housing. Parameters for all runs are listed in Table II. 

Group I was selected as the criterion for correlating the results in the several runs. These 
results are indicated in Table III and Figures 7 and 8. Table HI includes the slope A5^p/A( AT ) and 
the corresponding interaction loads for the respective runs. The longitudinal interaction force, 

X H , and interaction moment, X M , are shown respectively in Figures 7 and 8 for each run. The 
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results indicate that run 4 with Beryllium as the structural material yields the lowest deflection 
for a given transverse temperature gradient, AT. This run also lends itself to the highest in- 
teraction loads as might be expected. Comparing runs 1 and 2, the effect of zero thermal gradient 
across the experiment housing is to decrease the horizontal force, x H , by approximately 3% and 
increase the moment, X„, by about 13%. However, as displayed in run 3, a AT of 7°F greatly re- 
duces x M by about 80% while increasing x H by approximately 20%. It is of interest to observe that 
the deflection curve slope is not affected appreciably by changes in at in the range considered 

i.e., 0 - 13 °F. Considering runs 1 and 6, it is seen that increasing the experiment housing cross 
sectional moment of inertia by a factor of ten increases x M by approximately 37%, decreases x H by 
about 6%, and, most significantly, reduces the deflection slope by approximately 31%. Doubling the 
support tube length or, more meaningful, increasing the distance between experiment mounting 
points (run 5) yields a change in the deflection slope which is indicated to be directly proportional 
to this al. This relation is seen by dividing Equation (6) by at and recalling that x v equals zero. 
An interesting result in comparing runs 1 and 5 is that the interaction loads are not a function of 
L. This result can be shown by expanding Equation (4). This expansion also shows the vertical 
load, x v , to be zero as indicated in Appendix A. 


CONCLUDING REMARKS 

For the mathematical model of one experiment combined with the support tube under the in- 
dicated assumptions, the following points may be concluded from the results. 

1. It is possible to obtain the same misalignment or rotational deflection for different com- 
binations of t, t, and at. 

2. It is possible to maintain zero deflection for at in the range of ±7°F if the values of T, T and 
at are properly chosen. 

3. Changes in T and T are more significant to misalignment than changes in at. 

4. The rotational deflection is directly proportional to L. 

5. The interaction loads are independent of changes in L. 

References 

1. McMinn, S. J.: "Matrices for Structural Analysis," John Wiley and Sons, Inc., New York, 1962. 

2. "Design Study of AOSO RAC-762-5-1 Tech. Vol. I," Republic Aviation Corp., Long Island, New 
York, 1962. 

3. "Honeycomb Sandwich Design,” Brochure "E,” Hexcel Products Inc. 

4. Hoff, N. J., and Mautner, S. E., "Bending and Buckling of Sandwich Beams," Journal of the 
Aeronautical Sciences, Vol. 15, No. 12, pp. 707-720, December 1948. 


7 



APPENDIX A 


I. DEVELOPMENT OF INDIVIDUAL TERMS WITH EQUATION (2) 
A. Temperature induced displacement and rotation matrices 


r« 


■\ 







V T 


> 


& 
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T J 


(A.l) 



Sketch A-] 

Displacements due to a uniform 
temperature increase. 



Sketch A-2 

Displacements due to a linear vertical 
temperature gradient. 


From Sketches A-l and A-2, induced displacements at point (2) relative to point (1) due to a 
uniform temperature increase (T) and a linear vertical temperature gradient (AT) are shown re- 
spectively to be 


and 


hj = aTL , &J 
"t t 


8 / 
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; At = 


0 * , 


. At 
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2h • 



aAIL 
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•Assuming small deflections. 


(A. 2) 


(A. 3) 
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The induced displacements of point (P) relative to point (1) can be expressed as 


s u = s 


aATL 


s, 1 , * K) i ^ 


+ 


8 


„ aATL 2 
0 + 2h 


*e T = + » At = 0 


aATL 


Therefore, Equation (A.l) becomes 


a L 




!f ,\ ~ J aATL2 

{M - < 2h 


-aATL 


Going through a similar development, it can be shown that for point (P) relative to point (3), 

' (»-?)] 


(5,) 


aL 


aATL 2 

2h 

-aATL 


B. Influence coefficient matrix for loads applied at the neutral axis 


[f] 


f f f 

1 HH *VH 


f H V f VV f^V 


f ee 


(A. 4a) 
(A. 4b) 
(A. 4c) 

(A. 5) 


(A. 6) 
(A. 7) 
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Where f u is defined as the displacement in the j direction due to a load applied in the i direction 



Sketch A-3 

Unit loads applied at the 
neutral axis, point (2) . 

For unit loads applied at point (2), see Sketch (A-3), terms in Equation (A. 7) for displacements at 
point (2) become 


f HH AE f VH 


= 0 f 


6 H 


L 3 -L 2 

f HV " 0 f vv ’ 3EI f ev = 2EI 


f HS 0 


he ' 2EI t ee ~ El 


Therefore, for displacements at point (2) due to unit loads at point (2), Equation (A.7) becomes 


[f] 


JL_ 

AE 


0 

0 

L 3 

-L 2 

3EI 

2EI 

-L 2 

L 

3EI 

El 


(A. 8) 


Considering displacements at point (P) due to unit loads applied at point (2) as shown in Sketch 
(A-3), the following two changes occur in the individual terms of Equation (A.8): 
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VH 



L 2 h 
~ 4EI 


f 0H 



-Lh 
2 El 
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X v 

Sketch A-4 


Unit loads applied at point 
(P) off neutral axis. 


(A. 11) 
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With reference to Sketch (A-4), the indicated loads applied at point (P) reflect the following loads 
at point (2): 


= PX H 


x v = % 


X* = p x« - p x„ l 


Substituting these expressions into Equation (A. 11) and dropping the pre- superscript (P) yields the 
following load matrix at point (2) (neutral axis) due to off axis applied loads (point (P)): 






(A. 12) 


Similarly, for loads applied at point (P) (see Figure 2) referenced to the neutral axis, point (4) 


r 'n 


{x} = X x v 


(A. 13) 


X* + x H y 


U. IDENTITY FOR VERTICAL LOAD, X v 

In Part A of the results it was first shown that, for the assumed configuration, the vertical 
force, x v , is always zero. This result can be shown as follows: Writing Equation (4) in con- 
tracted form as 


{x} = [f +f]- 1 {S T + S T } 

and redefining the inverse of the combined influence coefficient matrix as 


(A. 14) 


[f +f]-‘ = [F ]' 1 = [k] 


(A. 15) 
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and the combined temperature induced displacement and rotation matrix as 


{S t + 8 t } " W 


Equation (4) can be also expressed in contracted form as 


(A. 16) 


M = [k] {d} 


(A. 17) 


[k] is known as the combined stiffness matrix. Expanding Equation (A. 17) yields 


Solving for x v 


X„ 


k ll 

k 12 

k 13 


S’ 'N 

d, 

Xy 

> = 

k 2l 

k 22 

k 23 

< 

d 2 

X, 

J 


_^31 

k 32 

k 33j 


d 3 

L J 


Xy kji dj + k 22 d 2 + k 23 d 3 


(A. 18) 


(A. 19) 


Where k 21 , k 22 and k 23 are the second row terms of [F] 1 in Equation (A. 15). 

The inverse of [F] is equal to the adjoint of [f] divided by the determinant of [F] . Since [f] 
is symmetrical, the transpose of [f] equals [F] . Further, the adjoint of [F] becomes just the 
cofactor of [F] . Therefore, Equation (A. 19) can be expressed in the form 


x v 





(A. 20) 


where are the cofactors of the f^ terms of Cf] , and I F | is the determinant of [f] . 
Expanding the second row cofactors of [F] 


C 2 1 ° ~(f 12 f 33 - Fia f 32 ) ; c 22 = +(f„ f 33 +f 13 f 31 ); 

^"23 ~ (Fll F 32 - Fj 2 F sl ) 

Substituting into Equation (A. 20) yields 

Xy “ |(" Fj 2 F 33 + F 13 F 32 ) dj + (F n F 33 - F 13 F 31 ) d 2 + (~ Fj j F 32 + F 12 F 31 ) d-jj- (A. 21) 
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Substituting terms of Equation (4) into Equation (A.21) 


Xy 




/a+A h 2 h 2 \ (1 +x\ l/h h\ 2 l faSL/AT At\ 

_\ AA + 41 + 41 / \ II / 4 \T l ) [_ 2 \h h /_ 

L/A+A hf hA L/E h\ 2 
2 V AA + 41 + «/ “ 8 \i - l)_ 


| aE[ jpL_^)l(A.22) 


Since C 21 d 1 = 0 and C 22 d 2 = -C 23 d 


3 ’- 


x v = 0 


(A. 23) 


for all cases. 
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APPENDIX B 


Fortran Listing of Programs I and n 


I. SYMBOL INTERPRETATION 

Fortran 

A 

AH 

AI 

AL 

DT 

DTB 

T 

TB 

EE 

ALPHA 

AIB 

AHB 

AB 

Z(ij) 

C(ij) 

D 


Problem 

A 

h 

I 

L 

AT 

AT 

T 

T 

E 

a 

T 

h 

A 

(f + fh, 

Cofactor of (f + f) Sj 
Determinant of [f + f ] 
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HI. LISTING OF PROGRAM R 


P 

C PKuGKhi i 1 I 

inmjihisM c t^b fiiqi I 

lu FukMAT (3ti4.4) 

l_i FuKPhT ( 1 HU »4A * 3HX ( H ) E_F_1 0# 4 / 3 X 9 bhK ( V ) = 1 U >4 /dX » 3 HX ( M ) = F ±0 . 4/ ) 

12 FukhmT ( LbC-lZ.4 } 

13 F uRhaT ( Inu ) 

14 FOki iAl ( 7F J.3,6 * iuX 94hC AbL I 4 } 

13 FUKhaT (^VX» 10Ftt>HT^ I NHUT ) 

16 FoRmAT ( lHU*30X*7HKt bUL T o ) 

17 FuknmT ( 1 nu »6 X » 1Ha» 1 1a » ihm 1 ux» Ittb? 1 uX-> 3h*LPhA * 6 X » ;>h 1 D ,, K »7X» 

33 MH cAK * 7X * 3HA dAK *yA»ihl »i0a* 1 HL • iOxmHL*-) 

lb F uKi-iaT ( Ifiu »7X » 7hULL 1 m 7»l2X»iHT»12X»3HT k » 7 X » 1 inu lLTm 1 p mk » 6_X * 

bbHULFHc: X * bX • 6HufcF RFL * 9X * feHu LF t kb ) 

1 kLhD INPUT ThPE 2»luf t~\ 9 i\Hi ll? ALPHA* A 1 3 • MtD* f U'« ^ 1 * AL* mL^ 
^ klJL uuTPoT TAPE 3,2^ “ “ “ 

UNITE. OUTPUT TaPL 3»17 

w K I T t 1 Pol lf-HL -5*l<i* ^ * 77 hT Ct« m(_PH/-W -\ i t.' * 77 I L 9 7c 9 M i y hL» A L £ 
wkiTt uuTHUT TAPE 3* lb 
iG^bb = 1 

UU 2 G J = 1 » 3 *2 

UT = j 
Nb 1 UH= 1 

O UO *0 K=1 »i^»4 

T=g*GSIGu 

Lb I GN= I 

9 L>U 33 L = 1 » 3 » 2 _ 

T b=L*Lb 1 GN 
nbIGN=l 

7 uu 23 rrrTTTTG 

uT b = i i*i ibl Oh 

TTTT1 = ( ho+^ ) / { aE*77 ) h77u:.. 2 / Jo)+h H * * 2 / ( 4 • -*- m I ) 

Z { 12 ) = -hL/4,* ( Ahlj/*u>-Ap/mI ) 

ZTT31 = • 3* ( M-ib/MTij-AH / * I ) ~~ 

Z ( 3 3 ) = ( m I b+A I ) / ( A I b* A I ) 

2(277 = "AL^^Z7"3.^rr3bJ 

Z(23) = -*L/2**Z(33) 

2- ( - i j = z < a k ) 

Z ( 3 1 ) = Z ( 1 3 ) 

/(3 ^j = z ( 2 3 ) 

C(li) = Z ( 22 ) *Z (33) — 2 ( 2 3 ) *Z ( 32 ) 

H I T T " 1 ^=TT 2 1 P TZ ( 33 TT7 ( 77 T ^T T TU 

U13) = Z(21 J*Z { 32 J-Z (22 )*Z ( 31 ) 

C < 2 1 ) = -Z ( i2 )*Z{ 33 )+Z( 13 )*Z<32) ‘ 

G ( 2 2 ) = Z(11)*Z(33)-Z(13)^Z(j>1) 

C ( 2 3 ) =-Z ( 1 1 )"*Z ( 32 ) +Z ( 1 2 ) *L ( 3 1 f 

C ( 3 1 ? = Z(l2)»Z(23)-Zil3)»Z(22) 

C ( 32 ) =-Z<lI)*Z(23)+Z(13)*Z(21) 

G ( 3 3 ) = ZUD*Z(22)- Z(l2)*Z(Zl) 

u = Z( 11 ) 22 }*Z ( 33 ) +Z ( 12 }*; (Z3 )*Z { 31 ) +7 :< 13T*k ( 3Z'T*Z Tz 1 ) 

3 -Z ( 13 )*a ( 22 )*Z ( 31 )-Z ( 12 )*2 { 21 )*Z ( 3 3 )-7 ( 1 1 > *Z ( 3 2 )*Z ( 2 3 ) 

Y = LL* ALPHA /U 

Z ( G 1 ) = (Tb-T-dHb+uJ )/2. )*Y 
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WRITE OUTPUT TAPE 3* il » ( F ( I ) 9 I = 1 , 3 ) 


utFHEX = 1 -aLPhm*uTd*AL/AHl-> + mL**^*F ( ^ ) / (?**ll*aIl>)--aL/ ( LL*m I b I 

( F ( 3 )+F ( 1 )*AHb/2. ) )*Z.06i6tb 

UtEkEL = -*LPhA*-uTo*AL,'*2 •06^6t5/Ahb 

utF LNL»= U l M i L X + u 1 1 RlL 


vv KITE O u l Pol TAPE o 9 1 A » \j I » T «► u I b ♦ n h F ri fc X * u L F K t L * U t F t IM u * 

& I C Abb „ 

^b COimT I fvUL 

I F ( UG I O INi ) o A » 3 U » 2 6 

lb i v iG I GPm = - 1 ” 

GO To 7 


34 I GaGE = I C A G t:. + 1 

3b CONTINUE 

I F ( L G I G i\ ) 4 b 9 3 0 * 3 6 

J 6 L G I G in = ~ 1 _ 

GO TO V 



APPENDIX C 


* 


Moment of Inertia Calculation for the Proposed Honeycomb Sandwhich Hexagonal Tube 



+y 


-z* 


Sketch C-l 

As indicated in Sketch C-l, the proposed support tube is composed of six honeycomb sandwich 
beams joined to form a symmetrical hexagonal cross section. Considering a typical cross section 
of the honeycomb sandwich structure as shown in Sketch C-2, the moment of inertia about its y-y 
axis can be expressed as 


b(t 3 -t/) 
” 12 


(C.l) 


where 

t Sandwich thickness 
t c Core thickness 
b Section width 
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Sketch C-2 


Since the cross section of the tube is symmetrical across both the Y and z axes, only one quadrant 
need be calculated and then multiplied by four to obtain the total cross sectional inertia of the tube 
about the Y-Y axis. 

Referring again to sketch C-l, the -Y/-Z quadrant is shown to be subdivided into two sections, 
A and B. Expressing the total inertia of the tube in terms of these two sections 

b 

x v r 4 2I Pw + ( Z i') 2A J (C-2) 

i - a 

where 

I yi Centrodial area inertia about the y-y axis of section i 

Z/ Normal distance from Y-Y axis to centrodial (y-y) axis of section i 

A a Effective cross sectional area of section i 

Expanding Equation (C.2) 

ly = 4[l ya + (z;) 2 A a + I yb + ( 2 ^') 2 A b ] (C.3) 



Sketch C-3 
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With reference to Equation (C.l) and Sketches (C-2) and (C-3), the centrodial area inertias for 
sections A and B can be expressed respectively as 


!y. = 12 (^^c 3 ) 

b. r -n 

(C.4) 

!,b = 12 (t c ) 3 ] 

From Sketch C-3, t' and t c ' can be expressed in terms of t and t c as 

(C.5) 

t ' = t sec 4 > 

and 

(C.6a) 

t c ' = t c sec <f> 

Substituting Equations (C.6) into Equation (C.5) 

(C.6b) 

b. sec 3 4 > 

«* - 12 (<■-•.’) 

(C.7) 

Assuming the effective area for bending to be only the face sheet cross section of the sandwich 
structure, the effective areas become (see Sketch C-l) 

A, = b. (t - t c ) 

(C.8a) 

\ * 2A. 

Substituting Equations (C.4), (C.7), (C.8a) and (C.8b) into Equation (C.3) 

Iy = 4 { b a [ 12 + ( Z a? (t-tc)] 

b. sec 3 <p l 

(C.8b) 

* 12 

(C.9) 


The following values for parameters in Equation (C.9) were used to evaluate I y . 
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t = 0.774 in. 
t c = 0.750 in. 
b a = 4.165 in. 
b b = 0.894 in. 
z a ' = 6.827 in. 

Z b ' = 3.413 in. 

<p = 30° 

Solving Equation (C.9) for these values yields an l y equal to 28.03 in 4 . 
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Table I 


Grouped Run Cases with Thermal Parameters. 


Group No 

Case No 

AT 

T T 

Group No 

Case No 

AT T T 

Group No 

Case No 

AT T T 

i 

1 

1 

1 1 

49 

49 

3 11 

97 

97 

5 1 1 


9 

1 

5 5 


57 

3 5 5 


105 

5 5 5 


28 

1 

-1 -1 


76 

3 -1 -1 


124 

5 -1 -1 


36 

1 

-5 -5 


84 

3 -5 -5 


132 

5 -5 -5 

2 

2 

1 

1 3 

50 

50 

3 1 3 

98 

98 

5 13 
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73 

3-11 


121 

5-11 


35 

B 



83 

3 -5 -3 


131 

5 -5 -3 

3 

3 

B 


51 

51 

3 1 5 

99 

99 

5 15 
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B 



74 

3-13 


122 

5-13 


34 

B 



82 

3 -5 -1 


130 

5 -5 -1 


42 

1 

-9 -5 


90 

3 -9-5 


138 

5 -9-5 

4 

4 

1 

1 -1 

52 

52 

3 1-1 

100 

100 

5 1-1 


8 

1 

5 3 


56 

3 5 3 


104 

5 5 3 


29 

i 

-1 -3 


77 

3 -1 -3 


125 

5 -1 -3 

5 

5 

1 

1 -3 

53 

53 

3 1 -3 

101 

101 

5 1-3 


7 

i 

5 1 


55 

3 5 1 


103 

5 5 1 


15 

i 

9 5 


63 

3 9 5 


111 

5 9 5 


30 

i 

-1 -5 


78 

3 -1 -5 


126 

5 -1 -5 

6 

6 

i 

1 -5 

54 

54 

3 1-5 

102 

102 

5 1-5 


10 

i 

5 -1 


58 

3 5-1 


106 

5 5-1 


14 

i 

9 3 


62 

3 9 3 


110 

5 9 3 

11 

11 

i 

5 -3 

59 

59 

3 5-3 

107 

107 

5 5-3 


13 

i 

9 1 


61 

3 9 1 


109 

5 9 1 


21 

1 

13 5 


69 

3 13 5 


117 

5 13 5 

12 

12 

1 

5 -5 

60 

60 

3 5-5 

108 

108 

5 5-5 


16 

i 

9 -1 


64 

3 9-1 


112 

5 9-1 


20 

i 

13 3 


68 

3 13 3 


116 

5 13 3 

17 

17 

1 

9 -3 

65 

65 

3 9-3 

113 

113 

5 9-3 


19 

i 

13 1 


67 

3 13 1 


115 

5 13 1 

18 

18 

i 

9 -5 

66 

66 

3 9-5 

114 

114 

5 9-5 


22 

i 

13 -1 


70 

3 13 -1 


118 

5 13 -1 

23 

23 

1 

13 -3 

71 

71 

3 13-3 


119 

5 13 -3 

24 

24 

i 

13 -5 

72 

72 

3 13 -5 


120 

5 13 -5 

27 

27 

i 

-1 5 

75 

75 

3-15 


123 

5-15 


31 

i 

-5 1 


79 

3-5 1 


127 

5-5 1 


41 

i 

-9 -3 


89 

3 -9 -3 


137 

5 -9-3 

32 

32 

i 

-5 3 

80 

80 

3-5 3 

128 

128 

5-5 3 


40 

i 

-9 -1 


88 

3 -9-1 


136 

5 -9 -1 


48 

i 

-13 -5 


96 

3 -13 -5 


144 

5 -13 -5 

33 

33 

i 

-5 5 

81 

81 

3-5 5 

129 

129 

5-5 5 


37 

i 

-9 1 


85 

3-9 1 


133 

5-9 1 


47 

i 

-13 -3 


95 

3 -13 -3 


143 

5 -13 -3 

38 

38 

i 

-9 3 

86 

86 

3-9 3 

134 

134 

5-9 3 


46 

i 

-13 -1 


94 

3 -13 -1 


142 

5 -13 -1 

39 

39 

i 

-9 5 

87 

87 

3-9 5 

135 

135 

5-9 5 


43 

i 

-13 1 


91 

3 -13 1 


139 

5 -13 1 

44 

44 

i 

-13 3 

92 

92 

3 -13 3 

140 

140 

5 -13 3 

45 

45 

i 

-13 5 

93 

93 

3 -13 5 

141 

141 

5 -13 5 
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Table II 


Parameters for Program II Runs. 


Run 

No. 

Material 

a 

(in/in/ °F) 

E 

psi 

Support Tube 
Length (in) 

L 

(in) 

L' 

(in) 

I 

(in 4 ) 

I 

(in 4 ) 

AT(°F) 

Range 

h 

(in) 

h 

(in) 

A 

(in 2 ) 

A 

(in 2 ) 

1 

Aluminum 

12.8 x 10- 6 

10 7 

60 

10.51 

16.17 

3.514 

28.03 

1 - 5 

7 

14 

0.5 

1.2 

2 

Aluminum 

12.8 x 10' 6 

10 7 

60 

10.51 

16.17 

3.514 

28.03 

0 

7 

14 

0.5 

1.2 

3 

Aluminum 

12.8 x 10" 6 

10 7 

60 

10.51 

16.17 

3.514 

28.03 

7 - 13 

7 

14 

0.5 

1.2 

4 

Beryllium 

7.39 x 1(T 6 

4.2X10 7 

60 

10.51 

16.17 

3.514 

28.03 

1 - 5 

7 

14 

0.5 

1.2 

5 

Aluminum 

12.8 x 10' 6 

10 7 

120 

48.00 

12.00 

3.514 

28.03 

1 -5 

7 

14 

0.5 

1.2 

6 

Aluminum 

12.8 X 1CT 6 

10 7 

60 

10.51 

16.17 

35.14 

28.03 

1 - 5 

7 

14 

0.5 

1.2 




Table III 


Deflection Curve Slopes and Interaction Loads for Program II Runs. 


Run No 

Case No 

AT(°F) 

AT(°F) 

Slope* 

(Arc Sec/°F) 

X„(Lb) 

X M (in-lb) 

1 

1 

+ 7 

+ 1 

-0.794 

- 97.18 

-369.21 

2 

1 

+ 7 

0 

-0.794 

- 93.82 

-418.47 

3 

1 

+ 7 

+ 7 

-0.794 

-117.34 

- 73.61 

4 

1 

+7 

+ 1 

-0.455 

-237.00 

-898.41 

5 

1 

+ 7 

+ 1 

-3.628 

- 97.18 

-369.21 

6 

1 

+ 7 

+ 1 

-0.547 

- 87.76 

i 

-507.37 


•Slope = AS 0 -/A(At) or DEFHEX/A(At) 
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